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Abstract 

A theoretical  investigation  has  been  carried  out  on  the  effect 
of  a drift  current  on  surface  polaritons  in  n-type  silicon.  The 
current  direction  is  taken  to  be  parallel  to  the  direction  of  propa- 
gation of  the  surface  polaritons.  From  the  dispersion  curves,  there 
is  evidence  that  an  interaction  takes  place  between  the  current  and 
the  polaritons  which  gives  rise  to  polariton  instability  or  ampli- 
fication for  certain  frequency  ranges.  These  instabilities  are 
believed  to  be  due  to  the  presence  of  the 
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I.  INTRODUCTION 


( 1— 5) 

Several  theoretical  investigations  ' have  been  made  of  the 
possibility  of  observing  surface  wave  instabilities,  e.g.  amplifying 
surface  waves,  in  a semiconductor  due  to  drifting  current  carriers. 
The  effect  of  a drift  current  on  surface  waves  has  also  been  observed 
experimentally v ’ . in  these  studies,  both  electrons  and  holes 

were  present  in  the  semiconductor.  In  the  present  investigation, 
however,  we  restrict  our  attention  to  only  one  type  of  current 
carrier.  We  have  theoretically  investigated  surface  polaritons  in 
n-type  Si  in  the  presence  of  drifting  electrons.  The  primary 
objective  of  this  investigation  was  to  determine  whether  or  not 
the  surface  polaritons  exhibit  instabilities  due  to  the  drifting 
electrons  when  only  electrons  and  not  holes  are  present. 

II.  SURFACE  POLAR ITON  DISPERSION  RELATION 


We  begin  by  using  non-local  Maxwell’s  equations  to  obtain 
equations  relating  the  field  and  dielectric  tensor  components. 

We  next  obtain  an  expression  for  the  dielectric  tensor,  which  is 
spatially  dispersive.  Imposing  the  so-called  dielectric  approxi- 
mation, partially  transformed  components  of  the  dielectric  tensor 
are  obtained.  Using  these  results,  the  above-mentioned  equations 
are  solved  for  the  field  components,  boundary  conditions  are 
applied,  and  the  surface  polar iton  dispersion  relation  is  obtained. 

Consider  a semiconductor  which  is  Infinite  in  the  x-  and 
y-directions  and  semi-infinite  in  the  z-direction.  In  seeking 
solutions  to  Maxwell’s  equations,  we  use  the  following  form  for 
the  arth  Cartesian  component  of  the  electric  field  E, 


m . 

wicmwxm  codes  1 

Dist 

C lldTu  IPLCMLI 

A 

3. 


Ea(x,t)  - Ea(z)ei(kxX  " u,t)  , (2.1) 

where  kx  is  the  wave  vector  and  uu  is  the  frequency.  Similar 
forms  obtain  for  the  displacement  field  D and  magnetic  field 
H.  In  this  paper  we  consider  only  surface  polaritons  with 
p-polarization,  i.e.,  the  electric  vector  lies  in  the  sagittal 
plane  defined  by  the  direction  of  propagation  and  the  normal  to 
the  surface.  Consequently,  Ey  - 0,  and  if  we  eliminate  the  magnetic 
field  from  Maxwell's  curl  equations,  the  latter  become 


"T5  Ex(z)  + ikx  ik  Ez(z)  ’ h Dx(z) 
dz  c 


(2.2a) 


lkx  Ex(z>  + kx  Ez<2>  - T Dz<z> 

c 


(2.2b) 


A functional  relationship  is  needed  between  D and  E to 
complete  the  set  of  equations.  For  the  crystal  (z  a 0),  we  can 
write 

00 

D*(z)  - E j*  dz'  €aa(kxtu|zz')E0(z')  . (2.3) 

e o 

Substituting  Eq.  (2.3)  into  Eqs.  (2.2),  we  obtain  the  following 
pair  of  coupled  integro-dif f erential  equations 


d2 


* 2 r 

Ex(z)  + ikx  di  Ez(z)  "^5  J dz/exx(kxaJ|zz')Ex(z') 

2 r 

+ ^ J dz/€xz(kxou|zz/)Ez(z/)  (2.4a) 


4. 


Hi  Ex<z>  + kx  Vz>  - ^ f dz'«zx<V’lzz'>Vz'> 

o 

2 * 

+ ^ f dz'€Zz(kxUJlZz/)EZ(z,)  (2.4b) 

c o 

which  determine  the  field  components  E (z)  and  E (z). 

X z 

In  order  to  solve  Eqs.  (2.4)  it  is  necessary  to  have  expressions 

for  the  components  of  the  dielectric  tensor.  We  proceed  by  assuming 

that  the  motion  of  an  average  carrier  in  the  semiconductor  is 

(8) 

governed  by  the  transport  equation v J , 

[V  + (V- V)V]  - - + | [E  + i V x H]  - vV  , (2.5) 

where  V and  v are  the  carrier  velocity  and  collision  frequency, 
respectively,  E and  S are  the  total  electric  and  magnetic  fields, 
respectively,  m and  N are  the  effective  mass  and  particle  density, 
respectively,  and  vP  is  the  electron  thermal  pressure  gradient. 

In  what  follows,  the  thermal  pressure  gradient  is  neglected,  i.e., 

VP  - 0 (this  is  equivalent  to  neglecting  the  thermal  motion  of 
the  current  carriers). 

We  linearize  Eq.  (2.5)  by  writing  V - VQ  + v,  E - EQ  + e, 

H - 3 + h,  * lere  V , H , and  E are  uniform  and  time-independent 

o o o o 

quantities  while  v,  h,  and  e are  position  and  time-dependent 
deviations  from  the  uniform  and  static  quantities.  These  linear- 
ized expressions  are  substituted  into  Eq.  (1.1)  and  terms  of  like 
order  on  both  sides  are  equated.  We  consider  the  situation  where 


/ 


5. 


H - 0.  Nothing  is  obtained  from  the  zero-order  equation,  but 
o 

for  the  first-order  equation  we  have 


v + (VQ  * 7)v  - J e + A (Vo  x h)  - vv  , (2.6) 


where  the  dot  indicates  time-differentiation. 

t 

Assuming  an  exponential  variation  exp[i(k*r  - o»t ) ] for  the 

first-order  quantities  v,  e,  and  h,  we  can  use  the  Maxwell  equation 

v x e - - — h to  eliminate  h from  Eq.  (2.10)  and  thus  obtain  the 
c 

result 

m[  u)v+  ivv-  (k* VQ)v]  - iq  [e Vo  x (kxe)]  . (2.7) 

Solving  this  equation  for  v gives  the  result 


where 


(2.8) 
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(2.9) 


Vox 


k V uj-(k  V +k  V ) 

z oy  v x ox  y oyJ 


and  a ■ muu  (uu  - k • VQ  + iv). 

The  particle  current  is  given  by  G - NV  - GQ  + g,  where 

- N V , g - N 7 + nV_ , N - N +n,  N is  the  mean  particle 
o o o o o o o 

density,  and  n is  the  deviation  of  the  particle  density  from  its 
mean  value. 

The  equation  of  continuity  can  be  written  as 

n + 7 • g . o , (2.10) 

which,  using  the  definition  of  g above,  becomes 


6. 


n + VQ  • vn  + Nq(v  . v)  - 0 . (2.11) 

The  quantity  n is  a function  of  position  and  time  which  we  express 
in  the  form  n(x,t)  - n (k ,ou )exp[  i (k*  r - cut ) ] . Substituting  into 
Eq.  (2.11)  we  find  that 

k*  v 

n - N - — . (2.12) 

° uu-k«VQ 

Substituting  this  result  into  the  equation  for  g above  gives  the 
result 

n.  vo(E.v)T 

g(k,uu)  - N v + — — . (2.13) 

L uu-k  • V J 
o 

The  electric  current  density  can  now  be  obtained  from  the 
relation 

3 - qg  ~ (2.14) 

In  component  form,  this  relation  can  be  expressed  as 


ja(k,uu)  - E aag  (k,uu)eg  (k,cu)  . (2.15) 

P 

Using  Eq.  (2.8),  we  find  that  the  nonlocal  conductivity  tensor 
a (k,uu)  is  given  by 


cag(k,uj) 


Noq 


(k,uu)  + — — 


oa_ 


(D 


-k*  V. 


E kYM^(k,uu) 


« 


(2.16) 


The  effective  nonlocal  dielectric  tensor  is  related  to  the  nonlocal 
conductivity  tensor  by 

t,aS<C-“,)  <2-17) 

where  €s  ia  the  background  dielectric  constant.  Note  that  the 
dielectric  tensor  is  spatially  dispersive  due  to  the  presence  of  terms 
containing  the  wave  vector  E. 


7. 


Before  proceeding  to  obtain  the  surface  polariton  dispersion 

(9) 

relation,  a simplification  will  be  made.  Maradudin  and  Mills'  ' 
have  applied  the  so-called  dielectric  approximation  to  the  type 
of  dielectric  tensor  considered  here.  The  simplification  consists 
of  assuming  that  this  tensor  depends  on  z and  z'  only  through 
their  difference,  as  is  the  case  for  an  infinitely  extended 
medium.  Surface  corrections  to  the  dielectric  tensor  are  thus 
neglected.  This  is  consistent  with  the  assumption  that  the  polariton 
decay  length  considered  here  is  numerically  large  compared  to  the 
distance  over  which  the  dielectric  constant  varies  near  the 
surface. 

For  each  of  the  dielectric  tensor  components  obtained  from 
Eq.  (2.17),  we  need  to  evaluate  the  partially  transformed  dielectric 
tensor  given  by 

r dk  _ ik  (z-z') 

€„B0VIzz'>  - J ^ 6ae(k,u,)e  , (2.18) 

— 00 

where  we  have  taken  k - V «■  V - 0.  The  results  are 

y oy  oz 

€ 

6xz(kx„|zz')  - -is4  - «IX(V>IZ,'>  <2-19b> 

OZ 

€zz(kxuulzz')  " P5)6(z_z/)  » (2.19c) 

where  6(z-z')  is  the  Dirac  delta  function,  and  where  we  have 
used  the  definitions 


3 f3 

xx(kx«lzz'>  - <«.-|f>5<z-z') 


as  (z-z ' 


(2.19a) 


m2aJ2(aJ-kxVox)(u)-kxVox  + iv)‘ 


8. 

(2.20a) 


e2 

^ 5 


, 2 2 2,  ...  . . 
€muju)(uj-kV  + iv) 

f."2*pVo2x'*-VM  + lv) 

e-"1DVOX 

‘“2<“-kxVox  + lv) 

<“>-kxVox> 

a,2 


(2.20b) 


(2.20c) 

(2 . 20d) 


(2. 20e) 


2 2 
and  ui  - 4rrN  q /m€ 

p O oo 

The  integrations  in  Eqs . (2.4)  can  now  be  carried  out,  with 
the  result  that  we  are  left  with  the  coupled  differential 
equations : 


-d2Ex(z) 

_ 


+ ik. 


dEz(z) 

dz 


«,2 


(‘.-I*) 


2 

Ex<z>  + if  “ 


d2Ex(z) 


+ i 


2 

(JU 


dEz(z) 

dz 


(2.21a) 


dE  (z)  „ 

lk-  nz — + kxEx(z>  • 1 H — Sr~  + *5  <«.-Es>Vz>  • <2-21b> 

c c 


dE  (z)  2 

uu  . xv  1 (JU 


It  is  assumed  that  for  the  crystal  (z  s 0), 


Ex(z)  - Ce~“Z,  Ez(z)  - De”az  , (2.22) 

where  C and  D are  constants.  Substituting  these  expressions  into 
the  coupled  differential  equations  gives  us  two  equations  relating 
C and  D.  Setting  the  determinant  of  the  coefficients  of  C and  D 
equal  to  zero  gives  us  an  expression  for  the  decay  constant  a: 


V 


9. 


(e-~  lf)[kx-^  «.-v] 

«. ‘if  ^)-  ^S-2V. 


(2.23) 


From  Eq.  (2.21b)  we  obtain  the  following  relation  between  C and 
D:  .2 

l 

c 


D - C 


c 


(2.24) 


Another  field  component  needed  is  D (z);  from  Eq.  (2.3)  and  using 

z 


Eqs.  (2.19)  we  obtain 


Dz(z)  - ip 


dE  (z) 

4 —35—  + <«--P5>Ez<z) 


(2.25) 


Consider  next  the  electric  fields  in  the  vacuum  (z  < 0) 

where  D (z)  = E (z).  Assuming  that 
a a 


at  z _ or  z 

Ex(z)  »Ae°  , Ez(z)  - B e ° 


(2.26) 


and  proceeding  as  for  the  case  of  the  crystal,  one  obtains  the 
following  results  for  the  decay  constant 

2 


a2  . k2  _ 

®o  x 2 
c 


(2.27) 


and  the  relation  between  A and  B 

ik 

B * 


(2.28) 


The  boundary  conditions  at  z - 0 can  now  be  applied. 

Continuity  of  Ex(z)  gives  the  result  C - A,  while  the  continuity 

of  D (z)  produces 
z 


10. 


- - — A i04ac  + i«C 

o 


S4> 

C 

kx-  A.  (e»-e5) 


(2.29) 


from  which  we  obtain  the  surface  polariton  dispersion  relation 


€ - 0c  - k 3 . 

ao  K 5 x 4 


(2.30) 


III.  RESULTS  AND  DISCUSSION 

Theoretical  dispersion  curves  have  been  obtained  for  surface 
polaritons  in  the  presence  of  drifting  current  carriers.  The 
effect  of  retardation  is  included,  but  the  effect  of  damping  due 
to  scattering  of  the  carriers  is  neglected.  The  latter  assumption 
is  reasonable  in  high  mobility  material. 

Calculations  were  made  using  values  of  parameters  appropriate 
for  n-type  Si.  The  value  of  the  background  dielectric  constant 
6^  was  taken  to  be  11.7.  The  density-of-states  effective  mass 

x * jjco  1/3 

m =»  (m^m^.  ) used  in  the  calculations  was  obtained  from  the 

values^*0)  m*  - 0.19m  and  m*  » 0.98  m . The  electron  mobility 

tot  o 

(111  3 2 

Pe  was  taken  to  bev  J 1.6  x 10  cm  /V  sec. 

In  our  calculations,  the  thermal  pressure  gradient  was 
neglected  (see  Eq.  (1.1).  This  is  valid  if  the  following  in- 
equality holds  ^). 

V„ 

iT-  < 1 » (3-D 

ox 

* * 2 1/3 

where  v - e/m  p is  the  damping  constant  and  V_**(ft/m  ) (3rr  N ) 

6 r o 

is  the  Fermi  velocity.  The  practical  upper  limit  for  the  election 

8 

drift  velocity  is  ~10  cm/sec,  which  corresponds  to  an  electric 


11. 


field  of  ~105  volts/cm,  obtained  from  the  equation  vox  = Eoxue  » 

where  the  value  of  ue  is  as  given  above.  For  n-type  Si,  it  has 

been  shown  that  the  mobility  of  the  electrons  doesn't  vary  with 

electric  field  up  to  10^  volts/cm^^.  Three  values  of  drift 

velocity  were  considered:  V =■  0.0001c,  0.001c,  and  0.01c. 

ox 

For  the  parameter  values  given  here  for  n-type  Si,  an 
upper  limit  for  the  carrier  concentration  Nq  can  be  calculated 

2 * l 

from  Eq.  (3.1).  'Taking  Vqx=»  0.0001c  and  uu  =2  oip  = (4rrNoe  /m  6^)  , 

18  “3 

we  obtain  the  result  N S 4 x 10  cm 

o 

In  making  the  calculations  of  the  dispersion  curves,  real 
values  of  uu/uup  were  assumed,  and  complex  values  of  the  wave 
vector  k = k-  + ikQ  and  the  decay  constant  a = or,  + ia„  were 

XX  Z 1 Z 

calculated  using  a digital  computer.  The  values  of  m/uu  were 

P 

taken  in  steps  of  0.01  or  smaller  for  the  range  0 < ia/ujp<1.0. 

( 13 ) 

The  values  of  kx  were  calculated  using  Muller's  method  , a 
process  of  successive  approximations  from  initial  guesses  (which 
need  not  be  externally  supplied). 

Figures  1 and  2 show  the  dispersion  curves  for  the  situation 
where  the  drift  velocity  V = 0.0001c.  Note  that  for  k,  > 0, 
there  is  one  branch,  which  for  yu  3 0.5  uup,  follows  the  light  line; 
then,  for  higher  frequencies  it  bends  over  and  moves  to  the  right 
of  the  light  line.  For  frequencies  yy  & 0.95  uup,  there  are  three 
branches  for  k^  < 0 in  Fig.  2 which  will  be  termed  negative 
branches  in  the  discussion  that  follows.  Two  of  the  three  negative 
branches  occur  because  of  the  presence  of  the  drift  velocity;  for  if 
V - 0,  there  is  only  one  negative  branch,  a mirror  image  of  the 

OX 

positive  branch  with  k^  > 0. 


12. 


At  the  lower  frequencies,  two  of  the  negative  branches 
(branches  a and  b in  Fig.  2)  have  large  values  for  |ck1/uu  | 
which  decrease  as  the  frequency  increases.  The  other  negative 
branch  (branch  c in  Fig.  2)  originates  at  the  origin  and  increases 
in  frequency  as  |k^|  increases.  Eventually,  a point  is  reached 
where  branch  c intersects  branch  a.  To  the  right  of  this  inter- 
section point,  there  exists  a continuation  of  branch  a indicated 
by  the  dash-dot  lines  in  Figs.  1 and  2,  which  consists  of  two 
degenerate  modes  whose  wave  vectors  form  a complex  conjugate  pair. 
As  the  frequency  increases  above  the  value  at  the  intersection 
point,  this  degenerate  mode  moves  to  the  right  and  finally  termi- 
nates at  the  vertical  axis  at  uj  = w . In  addition,  the  real  part 

p 

of  the  decay  constant  a becomes  zero  at  the  termination  point. 

The  remaining  negative  branch  in  Fig.  2 (branch  b)  also  moves 

in  closer  to  the  light  line  as  the  frequency  increases,  and  then 

terminates  at  a point  where  the  frequency  is  slightly  below  uj  . 

P 

At  this  end  point,  the  decay  constant  a becomes  pure  imaginary, 
i.e.,  we  no  longer  meet  the  condition  for  surface  polaritons.  We 
note  also  that  branch  b intersects  branch  a at  a large  value  of 

lck1/lDpl  * 

As  we  see  from  Figs.  1 and  2,  there  is  a frequency  range 
where  we  have  a degenerate  negative  branch  whose  wave  vectors 
form  a complex  conjugate  pair.  For  all  the  other  negative 
branches,  the  wave  vectors  are  real  except  in  the  immediate 
vicinity  of  the  point  where  branches  a and  b intersect.  A plot 
of  the  magnitude  of  the  imaginary  part  of  the  reduced  wave 


vector  |ck2/ujp|,  versus  frequency  is  shown  in  Fig.  3 for  the  ? “ 

degenerate  pairs.  The  curves  for  |ck„/uu  | start  at  k,  - 0 
for  tu  — 0.96  uj  , rapidly  increase  with  increasing  frequency 

r 

to  a maximum,  then  reverse  and  approach  k2  - 0 again  as  ou 
approaches  ujp.  In  addition,  k2  is  non-zero  over  a very  narrow 
frequency  range  about  uu  - 0.92  uj  where  branches  a and  b inter- 
sect  in  Fig.  2. 

Figure  4 shows  the  frequency  dependence  of  the  real  part  of 
the  reduced  decay  constant,  cof./uu  , versus  frequency  for  V -0.0001  c. 

* r OX 

The  branches  are  labelled  to  correspond  to  Fig.  2.  For  that  fre- 
quency range  where  the  wave  vector  is  complex  conjugate  (see  Fig.  2), 
the  decay  constant  is  also  complex  conjugate. 

Figure  5 shows  the  dispersion  curves  for  an  order  of  magnitude 
higher  drift  velocity,  i.e.,  V - 0.001  c,  in  a region  of  negative 

OA 

k^.  For  our  purposes,  the  region  of  interest  centers  on  the  point 
where  branches  a and  c intersect.  Again,  we  have  complex  conjugate 
values  of  k in  the  upper  of  the  intersecting  branches  to  the 
right  of  the  intersection  point.  There  is,  as  was  the  case  for 
V - 0.0001  c considered  previously,  a positive  branch,  which 
although  not  shown,  is  essentially  the  same  as  that  shown  in  Fig.  1. 

We  also  have  complex  values  of  k^  in  the  immediate  vicinity  of 
the  intersection  point  of  branches  a and  b. 

Figure  5 differs  from  Fig.  2 in  two  noticeable  ways:  in 
Fig.  5,  the  intersection  point  of  branches  a and  c occurs  at  a 

smaller  value  of  |k^|  than  in  Fig.  2,  and  the  uppermost  negative 
branch  (branch  b)  moves  to  a smaller  value  of  |k^j  before  termi- 
nating. As  before,  the  termination  point  occurs  when  the  decay 
constant  no  longer  satisfies  the  requirement  Rea  >0.  In  addition, 
the  point  of  intersection  of  branches  a and  b occurs  at  a much 
smaller  value  of  |kjJ  than  in  Fig.  2. 
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In  Fig.  6,  the  frequency  dependence  of  the  imaginary  part  of 
the  surface  polariton  wave  vector  for  V - 0.001  c is  presented 

It  shows  a behavior  somewhat  different  from  that  shown  in  Fig.  3, 

the  case  for  V - 0.0001c.  The  frequency  range  where  complex 
ox 

conjugate  wave  vectors  occur  is  larger  for  the  higher  drift 
velocity,  while  the  maximum  value  of  Ickg/cUp)  is  smaller. 

Figure  7 shows  dispersion  curves  for  the  highest  drift 
velocity  considered  here,  namely  V - 0.01c.  The  relative 

OX 

behavioral  changes  discussed  in  the  last  two  paragraphs  continue 

for  V - 0.01c,  except  for  one  obvious  difference:  branches  a and  b 
ox 

now  coincide  over  an  appreciable  frequency  range,  the  wave  vectors 
in  this  range  being  complex  conjugates  of  one  another.  As  the 
frequency  increases,  the  branches  separate,  the  lower  one  finally 
intersecting  branch  c in  the  same  manner  as  shown  in  Figs.  2 and 
5.  The  intersection  point  differs  from  those  of  Figs.  2 and  5 
in  that  it  occurs  at  a smaller  value  of  |k^|.  In  addition, 
branch  b in  Fig.  7 moves  to  a smaller  value  of  | | before 
terminating  than  is  the  case  in  Figs.  2 and  5.  The 
plot  of  |ck0/uj  I vs.  uj  is  shown  in  Fig.  8 for  V _ - 0.01c. 

Note  from  Figs.  3,  6,  and  8 that  the  frequency  range  where  complex 
conjugate  wave  vectors  occur  is  larger,  the  higher  the  drift 
velocity,  while  the  maximum  value  of  Ickg/ujpl  is  smaller. 

We  now  consider  in  detail  the  degenerate  modes  whose  wave 
vectors  form  a complex  conjugate  pair.  One  of  these  modes 
has  the  imaginary  part  of  the  wave  vector  negative,  so  that  the 
amplitude  increases  as  the  mode  propagates.  We  therefore  have 


15. 


an  apparent  instability.  This  instability  may  be  caused  by 
the  presence  of  the  surface.  Specifically,  the  important  factor 
leading  to  the  instability  may  be  the  large  electron  density 
gradient  that  occurs  in  the  vicinity  of  the  crystal  surface. 

So-called  convection  and  nonconvection  instabilities  occur 
in  plasmas ^15’ 16 ^ and  are  labeled  according  to  whether  or  not  the 
signs  of  the  group  velocities  (duj,/dkx)  of  the  interacting  modes 
are  the  same  prior  to  switching  on  the  interaction.  If  the  group 
velocities  have  the  same  sign,  then  the  interaction  can  result 
in  a convective  instability .which  may  be  interpreted  as  an 
amplifying  wave.  If,  however,  the  group  velocities  have  opposite 
signs,  the  interaction  can  result  in  a nonconvective  instability 
where  certain  real  wave  vectors  correspond  to  complex  frequencies 
or  in  an  evanescent  wave  where  certain  real  frequencies  correspond 
to  complex  wave  vectors. 

For  the  interaction  associated  with  the  intersection  of 
branches  a and  c,  we  shall  regard  the  unperturbed  modes  to  be 
the  surface  polaritons  and  the  electron  drift  current.  The 
latter  is  taken  to  have  a positive  drift  velocity  (group  velocity), 
whereas  the  former  has  a negative  group  velocity.  Since  the  two 
group  velocities  are  of  opposite  sign,  we  expect  that  there  may 
be  either  an  evanescent  wave  or  a nonconvective  instability. 

Our  calculations  show  that  tht  degenerate  branch  under  discussion 
occurs  in  a range  of  real  frequencies  with  complex  wave  vectors. 
Consequently,  the  growing  wave  of  this  degenerate  branch  is  an 


evanescent  wave. 
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We  now  turn  to  the  second  type  of  interaction  which  is 
associated  with  the  intersection  of  branches  a and  b.  In  this 
case  the  two  group  velocities  of  the  interacting  modes  have  the 
same  sign.  As  already  stated,  this  branch  involves  a range  of  real 
frequencies  with  complex  wave  vectors.  We  have  also  verified  that 
if  we  choose  the  wave  vector  real,  the  frequencies  for  the  de- 
generate  branch  are  complex.  The  growing  wave  in  this  second 
type  of  degenerate  branch  is  therefore  an  amplifying  wave  and 
corresponds  to  a convective  instability.  It  should  also  be 
emphasized  that  the  instabilities  reported  in  this  paper  occur 
in  the  presence  of  only  one  type  of  current  carrier  and  do  not 
require  two  types  of  current  carriers  such  as  electrons  and  holes. 

In  a recent  paper  £ 16 ^ Tajima  and  Ushioda  have  presented 
theoretical  results  on  surface  polaritons  in  polar  semi- 
conductors in  the  presence  of  a current.  Under  certain  cir- 
cumstances, they  find  Instabilities  associated  with  the  coupling 
of  surface  plasmons  with  surface  optical  phonons,  but  do  not 
report  instabilities  of  the  type  discussed  in  the  present 
paper.  Tajima  and  Ushioda  used  a different  method  of  handling 
spatial  dispersion  than  the  present  authors  and  did  not  in- 
vestigate negative  values  of  the  wave  vector,  so  it  is  not 
surprising  that  they  did  not  observe  the  instabilities  reported 


here. 
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FIGURE  CAPTIONS 

Fig.  1.  Dispersion  curves  for  surface  polaritons  in  n-type  Si 
(Vox/c  “°-0001>  • 

Fig.  2.  Dispersion  curves  for  surface  polaritons  in  n-type  Si 

(vox/c  “ 0.0001) . 

Fig.  3.  The  frequency  dependence  of  the  imaginary  part  of  the  wave 
vector  associated  with  surface  polaritons  in  n-type  Si 
(VQX/c  - 0.0001) 

Fig.  4.  The  frequency  dependence  of  the  real  part  of  the  decay 

constant  associated  with  surface  polaritons  in  n-type  Si 

<vox/c  - o.oooi) . 

Fig.  5.  Dispersion  curves  for  surface  polaritons  in  n-type  Si 
( VQJ^c  - 0.001)  . 

Fig.  6.  The  frequency  dependence  of  the  imaginary  part  of  the  wave 
vector  associated  with  surface  polaritons  in  n-type  Si 
( vox/c  - 0.001)  . 

Fig.  7.  Dispersion  curves  for  surface  polaritons  in  n-type  Si 

<Vox/c  - 0.01). 

Fig.  8.  The  frequency  dependence  of  the  imaginary  part  of  the  wave 
vector  associated  with  surface  polaritons  in  n-type  Si 

<vo*/o  - 001)- 
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